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Abstract 

The electromagnetically induced transparency (EIT) phenomenon in ear- 
lier universe is considered. We evaluated the elementary processes of the 
single scattering of photon on the hydrogen atom with the purpose of their 
use in the tasks of the radiation transfer theory. The additional function /, 
which depends on external conditions, is found. This function can be con- 
^ ,— i sidered as an adjustment of the optical depth that leads to the necessity of 

modernization of the escape probability Pij(rs) — > Pij {ts (1 + /))■ The 
numerical values of / for the different schemes of atom in three-level ap- 
proximation are given. It is found that the magnitude of function / could 
influence significantly on the formation of CMB in some partial cases. 

1 Introduction 

Electromagnetically induced transparency (EIT) is a well-known quantum inter- 
ference phenomenon that arises when coherent optical fields couple to the states of 
a quantum system flX) . Interference effects arise, because in quantum mechanics 
the probability amplitudes must be summed and squared to obtain the total tran- 
sition probability between the relevant quantum states. Interference effects can 
lead to profound modification of the optical and nonlinear optical properties of 
medium. In EIT the interference occurs between alternative transition pathways, 
driven by radiation fields. Interfering transition pathways can be induced by appli- 
cation of resonant laser fields to multilevel atomic systems. The most interesting 
example for this type of interference is the cancellation of absorption for a probe 
field, tuned in resonance to an atomic transition. EIT exhibits a laser-induced 
interference effect in the atom, which leads to an almost transparent medium. 
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In this paper we investigate the medium transparency constructed from hydro- 
gen atoms under the cosmological recombination conditions. The transparency 
of primordial medium was achieved in cosmological recombination era, when the 
universe was cooled and the electrons combined with the protons, forming hy- 
drogen atoms. After "recombination", photons were able to travel through the 
universe relatively unimpeded, and formed comprise the primordial background 
radiation that should been affect on generated hydrogen. 

The problem of line transfer in medium can be described quite accurately with 
the Sobolev approximation 0. Within the resonance the photon may be scattered 
in medium by the resonant line many times before it escapes from the matter. To 
account for resonant trapping in the Sobolev approximation, one should replace 
the normal spontaneous transition rate by the effective transition rate, modernized 
by the Sobolev escape probability. Theory of radiation transfer with the use of the 
Sobolev escape probability pij, where j is the upper level and i is the lower level, 
for the multilevel atoms have been described in J3]|. 

The "standard" recombination calculation refers to the calculation with the 
"effective three-level hydrogen atom" that is widely used and first derived in flU, 
0, where the "effective three-level hydrogen atom" means a hydrogen atom that 
includes the ground state, first excited state and continuum. In Q it was noted that 
effective three-level hydrogen atom, which would completely neglect all levels 
above n = 2, would be a hopeless approximation. Good accuracy is obtained 
by considering an n-level atom, where n is large enough. In practice, authors 
||3l found that a 300-level atom is more than adequate for the description of the 
microwave background. 

In this paper the other kind of three-level atomic system which includes only 
the bound states is considered. We describe this system as a zero approximation 
to demonstrate the possibility of the EIT in recombination era of earlier universe. 
The three-level hydrogen is chosen as a system, that can be adequately reduced 
to three levels when interaction with the pertinent electromagnetical fields are 
considered. The atomic selection rules require that two pairs of levels are dipole- 
coupled, while the transition between the third pair is dipole forbidden. This 
formulation allows to separate the resonant (one-photon) and non-resonant (two-, 
three-photon and etc) radiation escape processes based on different effects. As 
such a system, we consider a hydrogen atom with the ground and two excited 
states with the principal quantum number n < 3. 

In Fig. 1 the basic three-level systems are shown. All of the level schemes, 
discussed in this paper, can be reduced to one or other of these schemes. We 
label the levels |1 >, |2 >, and |3 >. The dipole- allowed transitions are between 
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states |1 > and |3 > and between states \2 > and |3 >. Classification of the 
schemes then depends upon the relative energies of the three states: Fig. 1 a) 
ladder (or cascade) scheme with E\ < E 3 < E 2 , b) lambda (A) scheme with 
Ei < E 2 < E 3 , and c) vee (V) scheme with E 3 < E\ and E 3 < E 2 . In a lambda 
or ladder scheme state 1 1 > is normally the ground state of the atom. The lambda 
scheme has a special importance due to the possible metastability of state |2 >. 

Recently, on the basis of O-JHl, the ladder scheme with |1 >= \ls >, |3 >= 
\2p > and |2 >= \3s > under cosmological conditions of cosmic microwave 
background (CMB) was considered [9]. It was found that the additional function 
/ is arised in optical depth rs- The maximal value of the / was found and is about 
1.5% in case of exact resonances (when the both detunings are equal to zero). In 
case of exact two-photon resonance, when frequencies of fields are close but differ 
slightly to the corresponding resonances and the total detuning is equal to zero, 
we found / ~ 0.95%. It was expected that these modifications should bring a 
strong impact to the determination of the key cosmological parameters ifTOl . 

In present paper we describe three-level systems: ladder system with fine 
structure splitting, vee scheme and A scheme for the hydrogen atom. We study the 
response of the three-level system on the external fields originating from photons 
emitted during the recombination period and evaluate the absorption coefficient, 
which can be applied to the radiation transfer theory. The physics of atom inter- 
action with photon fields can be understood through interfering pathways descrip- 
tion which corresponds to multi-photon process defined by power series expansion 
over fields amplitudes. 

2 Theory and basic equations 

New cosmological experiments on the CMB anisotropy investigations allow to 
know about processes in earlier universe in details and to understand the phys- 
ical processes, playing crucial role on extra large scales in universe. Detailed 
understanding of the recombination process is important for modeling the power 
spectrum of CMB anisotropics. CMB distribution can be obtained on the basis 
of the radiation transfer theory, which was strongly re-examined recently, see, 
for example, O, ffTTTl - [fTVll . We will rely on the work [|3l for the search of pos- 
sible application of the EIT phenomena to astrophysics. In accordance to [0 we 
omit the evaluation of the radiation transfer equations and restrict ourselves by the 
derivation of the additional terms to the "standard" absorption coefficient, which 
is included in the optical depth. 
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The Sobolev escape probability can be defined as 

Pij=Pij( T s)i (1) 

where j is the upper level and i is the lower level, r s is the Sobolev optical depth 
that can be presented, in the special case, in the form ts = t^t. Here k is the 
integrated line absorption coefficient, V is the velocity gradient which is given 
by the Hubble expansion rate H(z) and \ij is the central line wavelength. The 
monochromatic absorption coefficient or opacity is k = k(j>(uij), fij is the fre- 
quency for a given line transition and (j)(vij) is the normalized line profile. 

Absorption coefficient depends strongly on the external conditions and re- 
quires the particular consideration for the each case. In presence of an external 
field the opacity can be defined from the density matrix py as follows: 

k = 9g o I™{Pij}, ( 2 ) 

Lij 

where e is the vacuum permittivity, N is the number of atoms, dij is the dipole 
matrix element, fiy is the Rabi frequency and u}y is the frequency of the corre- 
sponding transition \j >— > \i >. 

Further we assume that during the recombination period in earlier universe 
the hydrogen atoms reach their ground states with emission of all the spectra of 
photons corresponding to the continuum-bound and bound-bound transitions in an 
atom. Consideration of the "atom-field" system can be restricted by the sponta- 
neous emission rates. Collisional excitation and ionization can be omitted because 
at the relevant temperatures and densities they are negligible 0. 

In this case the application of the density matrix formalism in three-level ap- 
proximation seems to be more simple and useful. The clear description of the 
density matrix theory and its applications can be found, for example, in |[T8l . lfT9l . 
As it was shown in [|9l the line profile can be separated out in the imaginary part 
of density matrix in approximation of weak fields and it was established 

Im{p tj } « 0(^) (1 + f(n 2 p ,n 2 c ,A p ,A c )) , (3) 

where <j>(fij) is the line profile which is multiplied by the corresponding Rabi 
frequency. The f(Qp, O^, A p , A c ) represents the dimensionless function which 
depends on the Q p , f2 c and detunings A p , A c (which are corresponded to the 
"probe" and "coupling" fields, respectively). 
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Thus for the three-level system the integrated line absorption coefficient can 
be presented in the form: 




4e 



[i + /(nj, j& a,, A,,)] . 



(4) 



The radiation transfer theory is based on the elementary processes of the elec- 
tron scattering on the atom at rest. In fact, dependence of the escape probability, 
Pij( T s), for the photon in the line wing due to the expanding universe can have 
the complex form. The Sobolev approximation works at a certain phase which is 
well known. In more complicated cases the diffusion approximation have to be 
applied [1 1J. However, we restrict ourselves by the task of taking into account the 
external fields in the single scattering of external photon. In more common case 
we can write 



where the optical depth Tg corresponds to the "standard" definition. Then we 
should replace py (t s ) by the p {j (rjj (l + /(fij, ftjj, A p , A c ))) 

Further we illustrate the importance and calculate the contribution of the func- 
tion /(Op, A p , A c ) for the different kind three-level systems, see Fig. 1. 

3 Derivation of the imaginary part of the density 
matrix element 

Consider the three-level system which is subjected to the two external coupling 
and probe fields. In order to describe the response of such system on a field that 
was formed during the recombination era, we apply the density matrix formalism. 
Taking in mind the continual spectrum of the cosmic microwave background, we 
can always pick up the resonant frequency even in such way that detuning will 
be vanish, and, thus, avoid consideration of the Doppler effect. Moreover, in 
equations we leave only spontaneous widths as a natural parameter of the spectral 
line due to the collisional excitation (relaxation) and ionization can be omitted at 
the relevant temperatures and densities 0). 

The time evolution of the density matrix elements is given by the master equa- 
tion (in atomic units): 




(5) 



d 



i < i\[H,p)\j > 




(6) 
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where \i > denotes the energy states of the atom, i = 1, 2, 3, H is the Hamiltonian 
of the system, H = H + H int , which can be presented in the rotating wave 
approximation like: 

H = wi|l >< 1| + u 2 \2 >< 2| + w 3 |3 >< 3| (7) 
H int = -- [tt p e-' lulpt \3 >< 1| + n c e- luJct \3 >< 2| + c.c] , 

where f2 p , f2 c are the Rabi frequencies of the probe and coupling field, respec- 
tively, the detunings from resonance are given by A p and A c , ui p and u c are the 
frequencies of the corresponding electric fields, cuj with i = 1,2,3 are the energies 
of the atomic states. The Rabi frequency is defined by f2 = d^E, where E is the 
electric field strength for the corresponding transition. Finally, detunings can be 
defined as A p = u; 3 — oj\ — u p and A c = u; 3 — w 2 — u c . 

The optical response of a medium to a resonant light is determined by the 
off-diagonal elements p 2 i, p^i and p 32 of the density matrix, for which p^- = p* { 
holds. In the master equation approach the effect of decoherence is taken into 
account by including damping terms, where it is assumed that the off-diagonal 
matrix elements decay with the respective rates T^. In the case of p 32 and p 3 i 
for the A scheme, this contains spontaneous emission from state 1 3 > to 1 1 > or 
1 2 >, respectively, with the total spontaneous emission rate T 3 = r 3i + T 32 . In 
the presence of damping, the time evolution of the coherences is then given by the 
Lindblad equations [|T9l : 

pij = -iujijpij - i [H int , p)ij - JijPij, i 7^ 3 (8) 
Pa = ~i [Hint,p}n + 2J ^ijPjj ~ ^2 ^JiP™- 

UJj >idi UJj <UJi 

Here r^- gives the rate per atom at which population decays from the level j to i, 
jij in absence of collisional processes can be represented as 

T^-Fi + r,-), (9) 

where Ti denotes the total decay width of the level \i >, i.e. 

r,= Yl (10) 

k(w k <u)i) 
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Thus for the three-level system, Fig. 1, the master equations for the density 
matrix element can be written in the form: 



Pl2 = -(712 + iUJu)pi2 - iQpP32 + i^cPl3 , 
p23 = -(723 + ^ 23 )P23 + iftc{PZi ~ P33) + iQpp21 , (1 1) 

P13 = -(7i3 + ^i 3 )pi 3 - in p (pu ~ P33) + i^cPn , 

where Q p corresponds to the probe field with 1 1 >— > |3 > transition 

and Q c corresponds to the coupling field with |2 >— > |3 > transition. 

Introducing new variables p 3 i = p3ie lbJpt , p-^2 = Pz2e tulct and p 2 i = p 2 ie l( ^ p+u!c ' t 
for the cascade scheme (rotating wave approximation) we recieve 

P12 = -(712 - i(\ + A c ))pia - itt p p32 + itt c p 13 , 

P23 = -(723 + «A c )p 2 3 + i^c(p22 ~ P33) + ^pP21 , (12) 

P13 = -(713 - «Ap)pi 3 + iQ p (pn - P33) + «^cPl2 , 

The Eqs. (PT2l) can be easily converted to the A and V schemes by the replacement 
A c — )• — A c and A p — > — A p , respectively. 

Full population of the ground state corresponds to the condition: 

Ph(0) = 1 (13) 
p 22 (0) = p 33 (0) = 0. 

Equations (fT2l) can be solved in the adiabatic regime or steady-state approxi- 
mation, i.e. p\j = that assumes approximation of the weak fields (the popula- 
tions of the states are changed weakly). Analyze of the solution of the Eqs. (fT"2l) 
for the different kind of schemes was given in a lot number of papers, see for ex- 
ample HI, flU-IEl, [fT8l . |fT9l . The response of the atoms to the weak probe field 
resonant to the 1 1 >— > |3 > transition is the quantity of interest for EIT. 

For the p 31 and p 32 (H scheme) matrix elements in the adiabatic regime we 
obtain: 

[(7i2 + ^)(732 + »A e ) +n 2 p ] 

p31 ~ ( 7l3 + iA p ) [( 7l2 + i<5)( 732 + iA c ) + fia] + ( 732 + iA c )W c ' 



P23 = 
P21 



( 7l3 + iA p )D? p + ( 723 + iA L ) [( 7l2 + z5)( 7l3 + iA p ) + ' 

( 723 + iA c )tt p tt c 

(713 + iA p )n 2 p + ( 723 + lA c ) [(712 + Z5)( 7 13 + iA p ) + 



where 5 = A p + A c . 

Application of the EIT phenomena to the astrophysics on the basis of ||6l-[[8l 
was done in [9]. For short we will not illustrate the series expansion over f2 p and 
Q c variables in assumption the smallness of the electric fields strength. Applying 
the same procedure like in O-JSl and separating out the Lorentz line profile the 
imaginary parts of the p 31 can be written in the form, see AH : 

« [1 + f(n p , n c , A v , A c )} . (15) 

A p + 7l3 

The analytical form of the function / is too cumbersome and further we will 
operate with numerical values. In principal, expression (fl5l) should be integrated 
over frequencies, but we omit these calculations (which we leave for the future 
works) and give only the estimates in assumption of the parametrical dependence 
of the function / on the detunings and Rabi frequencies. The matrix elements p 2 3 
and p 21 are not given due to their suppression by the factor Q p , as it should be. 
This suppression corresponds to the fact that we consider the full population of 
the lowest (ground) atomic level. 



4 Numerical evaluation 

The corresponding physical meaning of the function / can be found in [9| and 
can be understood via the optical depth Eq. ©. The "standard" way of the optical 
depth definition can be done through the integrated line absorption coefficient that 
is the characteristic of medium. The optical depth (see Eqs. (31), (33) in j3]|), 
originally, does not take into account the influence of the external field. How- 
ever, from optics it is well known that external fields should change the parame- 
ters of medium, including the spectral characteristics. The graphs 2 show which 
processes are included in function / that appeared as a result of external field's 
influence. Namely, the a) graph corresponds to the "standard" definition of the 
Sobolev escape probability. However, the existence of the second (coupling) field 
leads to the delay of the electron on the excited states, impedes the final recom- 
bination, due to the additional processes depicted on the graphs b)-d). Function 
/ is the dimensionless one and shows the branching ratio of the processes b)-d), 
i.e. the magnitude of this function reveals the contribution of the external fields in 
respect to the usual one-photon absorption process a). 

While the universe further expanded and cooled down the electrons and pro- 
tons tended towards formation of hydrogen atoms. The temperature at this era is 
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well known T e 4500 — 3000 K. After "recombination", the photons released 
were able to travel through the universe relatively undisturbed, and formed the 
primordial background radiation. Such photon environment (background) should 
have influence on the hydrogen atom. The field amplitudes for a circular polarized 
wave can be obtained from the (thermal-averaged) spectral energy density 



ce \E\ 2 _ 2hv? j Au ij 1 



47T C 2 ' 

gfe s T e — ]_ 



(16) 



where c is the speed of light, ks is the Boltzmann constant, h is Planck's constant 
and in further calculations we use T e = 3000 K. The right-hand side of the equa- 
tion above corresponds to the black-body distribution of the CMB, while left-hand 
side defines the (electrical) energy density. In order to avoid de-phasing problem 
we should choose Au^ ~ IV Thus, we can define the corresponding field ampli- 
tudes for the different H, A and V schemes by the fixing of the states |1 >, |2 > 
and 13 >. 



4.1 Ladder or cascade scheme S 

In BH the Ladder scheme with |1 >= [Is >, |2 >= \3s > and |3 >= \2p > was 
considered. It was established that the maximal value of the / is about 1.5% in 
case of exact resonances (when both detunings are equal to zero). In case of exact 
two-photon resonance, when frequencies of fields are close but differ slightly to 
the corresponding resonances and the total detuning is equal to zero, we obtained 
/ ~ 0.95%. In phis section we conider fine structure of the hydrogen levels. Due 
to the Lamb shift the \2p > state is lower then the \2s > and the corresponding 
levels we can set as: |1 >= \ls >, |3 >= \2p > and |2 >= \2s >, Fig. 1 a). 
The level widths are equal to T 2s = 8.229 s~\ T 2p = 6.26826 ■ 10 8 s" 1 and the 
transition frequencies are i/ 3i = 2.4674 • 10 15 and z/ 23 = 1.057911 ■ 10 9 in s" 1 . 
Hence for the spectral Ly a line z/ 3 i, setting Ai/ij = T 2p , we obtain 

E p w 0.000068802 V/m = 1.33799 ■ 10" 16 a.u., (17) 
E c w 0.0017496 V/m = 3.40242 • 10" 15 a.u., 
tt p w 9.96713- 10" 17 a.w., 
Q c w 1.02073 • 10~ 14 a.w. 

With the estimates Eqs. (fTTT ) the numerical results of the function /, see Eq. (fT5l ), 
for the different values of the detunings are presented in Table 1 . 
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Table 1. The numerical values of the function f(Q p , Q c , A p , A c ) for the differ- 
ent magnitudes of detunings are presented. In the first column the values of the 
f(Cl p , Cl c , A p , A c ) are listed, in the second and third columns the detunings are 
written. 



f (S^pi Ap, A c ) 


Ap s- 1 


A c s- 1 


-1.38066 ■ 10~ 4 








-1.38066 ■ 10" 4 


T 2s = 8.229 


-r 2s = -8.229 


8.28426 ■ 10~ 5 


r 2p = 6.26826 ■ 10 8 


-r 2p = -6.26826 • 10 8 


-8.12364 • lO" 6 


T 2s = 8.229 


T 2s = 8.229 


7.25115 ■ 10^ 13 


r 2p = 6.26826 ■ 10 8 


T 2s = 8.229 


3.62557- 10- 13 


r 2p = 6.26826 • 10 8 


r 2p = 6.26826 • 10 8 



From the Table 1 follows that the maximal value of the / is about 0.01% in 
case of the exact resonances (when the both detunings are equal to zero). In case 
of exact two-photon resonance, when frequencies of fields are close but differ 
slightly to the corresponding resonances and the total detuning is equal to zero, 
we found / m 0.0083%. Thus, in contrast to the results (3, the maximum value of 
the function / for the ladder scheme lies beyond required accuracy of calculations 
of the formation processes of the cosmic microwave background. We should note 
also that in case of Ai^- = T 2s the maximal value of function / is of the order of 
10~ 12 due to the smallness of the field amplitudes and we can conclude that EIT 
phenomena does not require the fine structure consideration for H scheme. 

4.2 Lambda scheme A 

As it was mentioned above the system of equations (fT2~l) for the A scheme, see Fig. 
lb), can be solved with the replacement A c — > — A c . Let us assume |1 >= \ls >, 

|2 >= \2s >, and |3 >= \3p > then u 21 = 2.4674 • 10 15 Hz, v 31 = 2.9243 • 10 15 
Hz, T 2 = T 2s = 8.229 s" 1 and T 3 = r 3p = 1.89802 • 10 8 s -1 . Assuming also 
Az/ 32 ps T 2s and Az/ 31 ps T 3p in Eq. (fl~6l , we obtain 

E p Pd= 2.45765 • 10~ 18 a.u., (18) 
E c Pd= 1.17805- 10~ 14 a.u., 
Q p ^ 7.33142 ■ 10~ 19 a.u., 
Q c Pd 2.08452 • 10" 14 a.w. 
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The numerical values of the function / Q c , A p , A c ) are given in Table 2. 

Table 2. The numerical values of the function f(£l p , f2 c , A p , A c ) for the differ- 
ent magnitudes of detunings are presented. In the first column the values of the 
f(£l p , fi c , A p , A c ) are listed, in the second and third columns the detunings are 
written. 



f (^p; Ap, A c ) 


Ap s- 1 


A c s- 1 


-0.00189829 








-0.00189829 


T 2s = 8.229 


T 2s = 8.229 


0.00113955 


r 3p = 1.89802 • 10 8 


r 3p = 1.89802- 10 8 


-0.000112064 


T 2s = 8.229 


-r 2s = -8.229 


3.29832 • 10~ n 


r 3p = 1.89802 • 10 8 


T 2s = 8.229 


1.64916 • 10" 11 


r 3p = 1.89802 • 10 8 


_r 3p = -1.89802- 10 8 



If we assume Au 32 ~ T 2s and Az/ 31 fa T 2s then Q p fa 1.52655 • 10~ 22 and 
tt c fa 2.08452 -10 -14 in atomic units. The amplitude of the field E c is the dominant 
and is equal to the field magnitude in previous case. Thus we will have the same 
results as in Table 2. 

In case Az/ 32 fa T 3p and Az/ 31 fa T 3p we get 

tt p fa 7.33142- 10~ 19 a.w., (19) 
Q c ^ 1.00111 ■ 10~ 10 a.w. 

Estimates Eqs. (fT9|) show that we can neglect by the terms which are higher 
order over Q p . However, the approximation Eq. (fl~5l ) is correct under the special 
conditions due to the violation of the relation tli/jij <C 1, namely f2 c /7 12 which 
is not a small parameter now. The solution p 3i takes the form: 

P3i f P Q2 ■ (20) 

7 13 + iA p + 

Taking into account the 7 12 = r 2s /2 ~ 10~ 16 a.w., one can find that the be- 
haviour of the matrix element p 3 i is not lorentzian for the vanishing detunings, 
Ap, A c — > or if Ap = A c (5 = 0). In opposite case when one of the detunings is 
large enough approximation Eq. (fl~4l) is kept. The values of function / for some 
detunings are listed in Table 3. 
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Table 3. The notations are the same as in Tables 1, 2. 



f (^p; ^c; Ap, A c ) 


Ap s- 1 


A c s- 1 


0.000761156 


r 3p = 1.89802 ■ 10 8 


T 2s = 8.229 


0.000380479 


r 3p = 1.89802 • 10 8 


-r 3p = -1.89802 • 10 8 


-9.04486 ■ 10~ 7 


T 2s = 8.229 


r 3p = 1.89802 ■ 10 8 



From the Tables 2 and 3 follows that the maximal value of the / is approx- 
imately 0.1% for the equal detunings, see Table 2, or is close to 0.1% for the 
estimates Eq. (|20l) in case of big A p (first and second lines in Table 3). 

Consider now fine structure of the A scheme. In this case we choose |1 >= 
| Is >, |2 >= \2s > and |3 >= \2p 3/2 >. Then u 31 « 2.46741 • 10 15 , v m « 
9.96903 ■ 10 9 , T 3 = r 2p = 6.26826 ■ 10 s s' 1 and T 2 = T 2s = 8.229 s~\ In 
order to estimate field amplitudes we apply again Eq. (fl~6l) with A 31 ~ T 2p and 
Az/ 32 ~ T 2s : 

O p s! 9.96632 -10~ 17 a.w., (21) 
ft c « 1.10204- 10~ 10 a.w. 

Numerical results of the evaluation of the function / are listed in Table 4 which 
shows that at such estimates Eqs. (|2TI) the contribution of the function / is negli- 
gible. When we define A 3i ~ T 2s and Az/ 32 ~ T 2s we got the analogous results 
due to the smallness of the field amplitudes. 

Table 4. The notations are the same as in Tables 1, 2. 



/ (ftp, Q c , Ap, A c ) 


Ap s- 1 


A c s- 1 


-1.60962 • 10~ 10 








-1.60962 • 10" 10 


T 2s = 8.229 


T 2s = 8.229 


9.65771 • 10~ n 


T 2 p = 6.26826 ■ 10 8 


T 2 p = 6.26826 ■ 10 8 


-9.46834 ■ 10" 12 


T 2s = 8.229 


_r 2s = -8.229 


8.45246 • 10~ 19 


r 2p = 6.26826 ■ 10 8 


T 2s = 8.229 



However, if we define A 31 ~ T 2p and Az/ 32 ~ T 2p we get 

ftp w 9.96632- 10~ 17 a.w., (22) 
ft c w 9.61828 • 10" 14 a.u. 
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Table 5. The notations are the same as in Tables 1, 2. 



f (^p; ^c? A p , A c ) 


A p s- 1 


A c s' 1 


-0.0121124 








-0.0121124 


T 2s = 8.229 


T 2s = 8.229 


0.00729051 


T 2p = 6.26826 ■ 10 8 


T 2p = 6.26826 ■ 10 8 


-0.000729015 


T 2s = 8.229 


-r 2s = -8.229 


6.43847- 10~ n 


T 2p = 6.26826 • 10 s 


T 2s = 8.229 


3.21924- 10" 11 


T 2p = 6.26826 ■ 10 8 


-r 2p = -6.26826 • 10 8 



and consequent numerical results which are illustrated in Table 5. 

Thus contribution of the / (Q p , f2 c , A p , A c ) reaches the value about 1%. This 
effect is within the accuracy of the modern computations of CMB and, therefore, 
should be taken into account. 

4.3 Vee scheme V 

In this section we consider vee V scheme depicted in graph Fig. lc). With the 
requirement the |1 >= \3p >, |2 >= \2p > and |3 >= \ls >, we have z/ 23 = 
2.4674 ■ 10 15 Hz, u 13 = 2.9243 ■ 10 15 Hz, T 2 = T 2p = 6.26826 ■ 10 8 s" 1 and 
r\ = T 3p = 1.89802 • 10 8 s _1 . Assuming Az/ 23 w T 2p and Au 13 w T 3p in Eq. 
(TT6l) . we obtain 

tt p w 7.33142- 10~ 19 a.w., (23) 
Q c w 9.96713 ■ 10~ 17 a.w., 

and numerical values of the function / (fl p , fi c , A p , A c ) are given in Table 6. 

Finally, lets consider the fine structure for the V scheme, Fig. 1 c). We choose 
|1 >= \2p 3 / 2 >, |2 >= \2pi/ 2 > and |3 >= \ls >, Then we have z/ 23 = 
2.4674 • 10 15 Hz, u 13 = 2.46741 • 10 15 Hz, T 2 = T 2pi/2 = 6.26826 • 10 s s^ 1 and 
r\ = r 2ps/2 = 6.26826 • 10 8 s _1 . Assuming A^ 3 w T 2p and Az/ 23 w T 2p in Eq. 
(fT6l) . we again obtain the vanishing values for the function / (O p , f2 c , A p , A c ) of 
the order 10~ 16 . Thus from Table 6 follows that such processes are unimportant 
and we can fully neglect by the consideration of V scheme. 
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Table 6. The notations are the same as in Tables 1, 2. 



/ (fip, n c , a p , a c ) 


A p s- 1 


A c s' 1 


-5.70837- HT 16 








4.8248 • 10" 16 


T 3p = 1.89802 • 10 8 


r 3p = 1.89802- 10 8 


4.52944 ■ 10" 16 


r 3p = 1.89802 ■ 10 8 


_r 3p = -1.89802- 10 8 


4.14603 • 10~ 16 


T 2p = 6.26826 ■ 10 8 


T 2p = 6.26826 ■ 10 8 


-1.14167- 10- 1S 


T 3p = 1.89802 • 10 8 


T 2p = 6.26826 • 10 8 


1.76641 ■ 10" 17 


T 2p = 6.26826 ■ 10 8 


-r 2p = -6.26826 • 10 8 



4.4 Full settling in the 2 s state 

The 2s <H- Is two-photon decay process is able to substantially control the dynam- 
ics of cosmological hydrogen recombination @J, 0, allowing about 57% of all 
hydrogen atoms in the Universe to recombine through this channel. Taking into 
account this fact, we demonstrate briefly the influence of the EIT phenomena on 
the schemes Fig. 1. With the initial condition p 2s (0) = 1, but others are equal to 
zero, and estimates Eqs. (PTTl) for the cascade scheme we receive numerical values 
listed in Table 7. 



Table 7. The notations are the same as in Tables 1, 2. 



/ (fip, n c , a p , a c ) 


A p s- 1 


A c s- 1 


-0.000276132 








-0.000276132 


T 2s = 8.229 


T 2s = -8.229 


-0.0000552386 


T 2p = 6.26826 • 10 8 


-r 2p = -6.26826 • 10 8 


0.00003314 


-r 2p = -6.26826 • 10 8 


T 2p - T 2s = 6.26826 • 10 8 


-0.0000162473 


T 2s = 8.229 


T 2s = 8.229 


7.25115 • 10~ 13 


T 2p = 6.26826 • 10 8 


T 2p = 6.26826 • 10 8 



For the A scheme for the matrix element p 32 we get the negligible values of 
the function / (Q p , fl c , A p , A c ), where the maximal value of the order 10~ 12 is 
reached at the zero detunings and Eqs. (fl~8l) . Taking into account the fine structure 
with Eqs. (|2~2|) the maximal value of the function / was found on the level 10 -8 . 
We should stress that for the case of the full population of the 2s state the matrix 
element p 31 is suppressed by the magnitude of the factor Vt p Vt 2 c and, therefore, 
negligible. 
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5 Conclusions 



In the last decade the research of the CMB anisotropy was carried out by many 
authors, see, for example, IfTTi - ffTTll where the hydrogen atom with the set of states 
about 300 levels was described. Such number of states is required in order to 
receive the accuracy about 0.1%. After recombination photons were able to travel 
through universe and formed the CMB. Generated by this way field should affect 
the primordial atoms. In this case we are confronted with the consideration of 
the phenomena of electromagnetically induced transparancy which can lead to 
profound modification of the optical and nonlinear optical properties of medium. 
In this work we limited ourselves to the consideration of the EIT phenomena in 
respect to amendment of the medium characteristic, such as absorption coefficient, 
assuming its application to the radiation transfer theory that we leave for future 
works. 

We employed the quantum optical evaluation of the integrated line absorp- 
tion coefficient that assumes application of photon beams (laser). In astrophysics 
the photon beam diffusion is accounted for via the optical depth. Optical depth 
expresses the quantity of light removed from a beam by scattering or absorption 
during its path through a medium. The optical depth definition (see Eqs. (31), 
(33) in (31), originally, does not take into account the influence of the external 
field. However, from optics it is well known that external fields should change 
the parameters of medium, including the spectral characteristics. Accordingly, we 
derived the absorption coefficient from the imaginary part of the density matrix el- 
ement and the magnitudes of fields were deduced from the CMB distribution. We 
considered the different schemes in three-level approximation depicted on Fig. 1 
for the hydrogen atom. The corresponding physical meaning of the function / Eq. 
([TBI can be found, for example, in flU, see also O-JSl. 

The radiation transfer theory is based on the elementary processes of the elec- 
tron scattering on the atom at rest. In fact, dependence of the escape probability, 
Pij(rs), for the photon in the line wing due to the expanding universe can have 
the complex form. The Sobolev approximation works at a certain phase which 
is well known. In more complicated cases the diffusion approximation have to 
be applied ifTTll . However, we restricted ourselves the task of taking into account 
the external fields in the single scattering of external photon. We found there 
is the additional function / which depends on external conditions (field ampli- 
tudes and detunings from resonant frequencies). Introduction of the function / 
leads to the necessity of modernization of optical depth Eq. © and, therefore, 
Pij( T s) — > Pij (t s (1 + / (fi p , f2 c , Ap, A c ))). In principal, the expression Eq. (fT5l) 
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should be integrated over frequencies. However, in this paper the problem was 
limited to the finding the maximal values of /. The values of function / for dif- 
ferent schemes of the hydrogen atom in three-level approximation are listed in 
Tables 1-7. 

The cascade scheme, Fig. 1 a), was considered in [0 and it was found that 
the maximal value of the / amounts to about 1.5% in case of exact resonances 
(when both detunings are equal to zero). In this paper, in assumption of the full 
population of the ground state, we found that the maximal value of the function / 
is acheived for the A scheme with the fine structure of levels, when the states are 
defined as |1 >= \ls >, |2 >= \2s >, |3 >= |2p 3 / 2 > and field amplitudes are 
defined via Ai/ 31 = Az/ 32 = T 2s - The maximal value of / is of the order 1.2%, see 
Table 5. The other significant example of the EIT phenomena is the case when 
|1 >= | Is >, |2 >= \2s >, |3 >= \3p > and Az/ 32 = T 2s , Az/ 31 = T 3p with 
the estimates Eqs. (fT8l) . Under these conditions it is obtained that the order of the 
function / is about 0.2%, see Table 2. From Table 3 follows that the function / is 
close to the 0.1% for the vanishing detunings and estimates Eqs. (fl9l ). For all other 
cases we got the smaller values. However, on the basis of the work [9] and these 
results we can conclude that the modifications of the integrated line absorption 
coefficient, could influence significantly on the formation of CMB, Eq. Q. 

Moreover, we examined the particular variants of the EIT phenomena mani- 
fistation with the full population of the 2s state which is caused by the fact that 
57% of hydrogen atoms recombine through the 2s state @|. However, it is was es- 
tablished that in this case the EIT phenomenon is not essential at the modern level 
of the experimental accuracy. For example, consideration of the cascade scheme 
leads to the contribution about 0.03% (Table 7) and is more negligible in all other 
cases. 

The problem of dephasing appears when defining the magnitude of the field. 
The effect of dephasing leads to the spectral line broadening. To prevent dephas- 
ing phenomena we confined the definition of the field by a narrow strip. In this 
case the width of the corresponding line appears as a natural parameter. In our cal- 
culations we used the following relation Ai/y = A p (A c ). In addition we should 
note also the exponential behavior over temperature for the field amplitudes. With 
increasing of the temperature T e larger values for the amplitudes could be ob- 
tained, see Eq. (fT6l) . Accordingly, the contribution of the EIT effect can be more 
significant for higher temperatures. 
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I2> 



a) b) c) 



Figure 1: Basic three-level schemes: a) ladder or cascade (S) scheme with E\ < 
E 3 < E 2 , b) lambda (A) scheme with E\ < E 2 < E 3 , c) vee (V) scheme with 
E 3 < E x and E 3 < E 2 . 
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c) 
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d) 



Figure 2: The transition processes, occurring in the three-level ladder system and 
corresponding to the different term of the imaginary part of Eqs. (fT5l) . are de- 
picted: part a) of the figure represents the one-photon absorption processes which 
corresponds to the "standard" definition of the optical depth; part b), c) and d) 
represent the terms of the next orders in fl p and f2 c appearing due to the presence 
of the external fields and derived via series expansion of the matrix elements Eqs. 
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